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00 (54) Title: METHOD FOR PREDICTING GRAIN SIZE DISTRIBUTION FROM RESERVOIR THICKNESS 

(57) Abstract: The present invention is a method for predicting the grain size distribution at a designated location within a water- 
lain sedimentary deposit. Initially, the vertical thickness of the sedimentary deposit at the designated location must be determined, 
as well as the vertical thickness and grain size distribution at a second location different from the designated location. Second, a 
Q distance parameter corresponding to the two locations must be determined. Finally the distance parameter is used, along with the 
initially determined vertical thickness at both locations and the grain size distribution at the second location to calculate the grain 
size distribution at the designated location. 
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METHOD FOR PREDICTING GRAIN SIZE DISTRIBUTION 

FROM RESERVOIR THICKNESS 



5 FIELD OF THE INVENTION 

[0001] This invention relates generally to the field of geophysical prospecting 
and, more particularly, to seismic and sedimentary data interpretation. Specifically, 
the invention is a method for predicting the grain size distribution at any point within 
10 a water-lain sedimentary deposit based on the thickness of the deposit at that point 
and a measurement of grain size distribution and thickness at a single other point in 
the deposit. 

BACKGROUND OF THE INVENTION 

[0002] In the oil and gas industry, seismic prospecting techniques commonly are 
15 used to aid in the search for and evaluation of subterranean hydrocarbon deposits. A 

■ 

seismic prospecting operation consists of three separate stages: data acquisition, data 
processing, and data interpretation, and success of the operation depends on 
satisfactory completion of all three stages. 

[0003] In the data acquisition stage, a seismic source is used to generate an 
20 acoustic impulse known as a "seismic signal" that propagates into the earth and is at 
least partially reflected by subsurface seismic reflectors (such as, interfaces between 
underground formations having different acoustic impedances). The reflected signals 
(known as "seismic reflections") are detected and recorded by an array of seismic 
receivers located at or near the surface of the earth, in an overlying body of water, or 
25 at known depths in boreholes. The seismic energy recorded by each seismic receiver 
is known as a "seismic data trace." 

[0004] During the data processing stage, the raw seismic data traces recorded in 
the data acquisition stage are refined and enhanced using a variety of procedures that 
depend on the nature of the geologic structure being investigated and on the 
30 characteristics of the raw data traces themselves. In general, the purpose of the data 
processing stage is to produce an image of the subsurface from the recorded seismic 
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data for use during the data interpretation stage. The image is developed using 
theoretical and empirical models of the manner in which the seismic signals are 
transmitted into the earth, attenuated by subsurface strata, and reflected from geologic 
structures. 

5 [0005] The purpose of the data interpretation stage is to determine information 
about the subsurface geology of the earth from the processed seismic data. The 
results of the data interpretation stage may be used to determine the general geologic 
structure of a subsurface region, or to locate potential hydrocarbon reservoirs, or to 
guide the development of an already discovered reservoir. 

10 [0006] At present, results from the data interpretation stage are generally limited 
to the size and shape of large subsurface structures. The descriptions may, for 
example, provide an indication of the total volume of hydrocarbons which might be 
retained in such structures. However, seismic data does not provide direct 
information about the grain size distribution of reservoir sediments, such as sands. 

15 Nor does present technology allow the analyst to be able to determine preferred 
locations within a structure for drilling wells to efficiently produce the contained 
hydrocarbons, except in a very broad sense. When an exploration well has been 
drilled, grain size samples can be taken, but present technology does not allow an 
analyst to accurately characterize grain size in locations other than in the most 

20 immediate region of any such well. 

[0007] As will be understood from this summary of the data interpretation stage 
of the seismic prospecting operation, it would be desirable to be able to predict the 
grain size distribution within a subsurface structure of interest using only seismic data 
and a minimal number of exploration wells, preferably just one exploration well. 

25 Such a capability would facilitate estimation of the hydrocarbon volume in place and 
potential production rates early in the hydrocarbon exploration and development 
process. As will be understood to those skilled in the art, hydrocarbon volumes and 
production rates depend on a variety of factors, including the grain size distribution of 
the sedimentary deposit that makes up the reservoir in which the hydrocarbons are 

30 found. The present invention is a method that may be used for predicting the grain 
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size distribution at one or more locations including every location within a 
sedimentary deposit without having to drill multiple wells to obtain actual samples of 
the reservoir sediments. 

SUMMARY OF THE INVENTION 

5 [0008] The present invention is a method for predicting the grain size distribution 
at a designated location within a water-lain sedimentary deposit. The vertical 
thickness of the sedimentary deposit at the designated location is determined, as well 
as the vertical thickness and grain size distribution at a second location different from 
the designated location. A distance parameter corresponding to the two locations is 

10 determined. The distance parameter is used, along with the initially determined 
vertical thicknesses and the grain size distribution, to calculate the grain size 
distribution at the designated location. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0009] Figure 1 depicts the plan view boundaries and the velocity components of 
15 a fluid flow which may create a deposit. 

[0010] Figure 2 depicts an elevation of the fluid flow of figure 1 and the elevation 
profile of the deposit that is created. 

[0011] Figure 3 depicts a flow chart of the calculation process for an embodiment 
of the method of the present invention. 

20 [0012] Figure 4 depicts a plan view of a laboratory-generated deposit used in an 
example application of the method of the present invention. 

[0013] Figure 5 compares measured with predicted results for the example of 
figured 

[0014] Figure 6 depicts a plan view of the Wax Lake delta area in Louisiana. 

25 [0015] Figure 7 compares measured with predicted results for the Wax Lake delta 
example of Figure 6. 
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DETAILED DESCRIPTION 

[0016] In the following detailed description, the invention will be described in 
connection with its preferred embodiments. However, to the extent that the following 
description is specific to a particular embodiment or a particular use of the invention, 
5 this is intended to be illustrative only. Accordingly, the invention is not limited to the 
specific embodiments described below, but rather, the invention includes all 
alternatives, modifications, and equivalents falling within the true scope of the 
appended claims. 

[0017] The present invention is a method that may be used for predicting the grain 
10 size distribution at one or more locations including every location within a 
sedimentary deposit using seismic data and without having to drill multiple wells to 
obtain actual samples of the reservoir sediments. Specifically, the inventive method 
can be applied to find the grain size distribution at any point within a water-lain 
sedimentary deposit, such as sand, based on the thickness of the deposit at that point 
15 and a measurement of grain size distribution and thickness at a single other point in 
the deposit For example, the thickness of the deposit can be determined from seismic 
data or from an exploration well and the measurement of the grain size distribution 
can be determined from data from an exploration well. 

[0018] The present invention improves on prior art technology in several 
20 important ways. First, the present invention relies on an explicit treatment of the 
physics of the sediment transport mechanisms that generate subsurface reservoir 
deposits. That explicit reliance on sediment transport physics decreases the need to 
rely on empirical correlations to specify deposit properties, an approach which has 
been heavily relied on in the prior art. (e.g., Fournier, F., "A statistical methodology 
25 for deriving reservoir properties from seismic data," Geophysics Vol. 60, 1995, pp. 
1437-1450. Haldorsen, H. H. and E. Damsleth, "Stochastic Modeling." J. Petroleum 
Technology. April 1990, pp. 404-412.) Second, the present invention can be applied 
to both individual deposits and to the thickness and vertically-averaged grain size 
distributions of larger stratigraphic units. Prior art methods have generally not been 
30 extendable to larger stratigraphic units as a result of the multi-component nature of 
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such units, (e.g., Scheidegger, A. E. and P. E. Potter, "Downcuirent decline of grain 
size and thickness of single tufbidite beds: A semi-quantitative analysis." 
Sedimentology Vol. 17, 1971. pp. 41-49. Syvitski, J. P. M., and J. M. Alcott, 
"Grain2: Predictions of particle size seaward of river mouths." Computers & 
5 Geosciences, Vol. 19, 1993, pp. 399-446.) Third, the present invention may be used 
to predict a complete grain size distribution. Prior art approaches have generally been 
limited to prediction of a single parameter of the distribution, such as the median or 
the largest expected grain size (e.g., Sadler, P. M., "Bed-thickness and grain size of 
turbidites." Sedimentology 'Vol. 29, 1982, pp. 37-51.). 

10 



Introduction 

[0019J Table 1 and Figures 1 and 2 define certain parameters, which will assist in 
fully appreciating the present inventive method. The method is built on the 

15 assumption that the sediment, both in the flow and in the deposit, may be 
characterized by bins of a common nominal grain diameter. Each such bin may, for 
example, be characterized by a nominal diameter du and include grains having an 
actual diameter with typical ranges from 84% to 119% of d ( . Other bin range 
definitions are within the scope of the present invention. The sediment grains are 

20 deposited from a sediment laden fluid flow immediately above the deposit, where that 
fluid flow has a height h(xy) and a vector velocity Vfoy). The scalar components of 
the vector velocity V are u(xy) and v(xy). Velocity V(xj/) is assumed to be constant 
with depth within the fluid flow. The quantity of sediment in the fluid flow may be 
defined by a grain volume fraction Cfey) for each grain size bin. The grain volume 

25 fraction may be assumed to be constant with depth within the fluid flow. 

[0020] Figure 1 depicts a plan view of a fluid flow 10, with boundaries 12 and 14, 
in an (pc,y) coordinate system. At any given point (x,y) within flow 10, vector velocity 
V will have scalar velocity components u(?cy) and v(xy). As depicted in elevation 
view Figure 2, the height of fluid flow 10 at any point is given by h(xy\ and the 
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thickness of the deposit at any given point is given by z(x,y). Deposit thickness z(x,y) 
results from grains in fluid flow 10 depositing from fluid flow 10 over a time duration 
To. 



5 Table 1 



Height of Fluid Flow 


Kxy) 


Velocity (vector) of Fluid Flow 


V{x,y) 


Components (scalar) of Fluid Flow velocity vector 


u(xy), v(* j;) 


XT/vm-tTml niamAtpr o*F Qp^imf*Trf" {"rffllYIQ fvf it'll OtJMTI 

Size bin 


d, 






Vnlump FrartinTi rvf OraTiw of the ith CrTain Size bin in 
the Fluid Flow. The sum of this term over all grain 
size bins gives the total fraction of the flow which is 
composed of sediment 




Thickness of Deposit resulting from deposition 




Grain Size Distribution in the deposit Gives the 
fraction of total grain volume contained in the ith grain 
size bin. The sum of this term over all grain size bins 
is one. 




Grain diameter at the Nth size percentile of the grain 
size distribution in the deposit at location x y. 




Still-water settling velocity of sediment grains of 
diameter d{ 




Effective Settling Velocity of Sediment Grains of 
diameter di 


v,(dd 


Sediment Density (typically 2650 kg/m 4 ) 


Pi 


Fluid Density (typically 1000 kg/m') 


Pw 


Gravitational Constant (9.8 mis*) 


g 


Viscosity Coefficient of Fluid 

(for water, typically 0.000881 Ns/m 2 ) 


M 


Time Duration over which deposition occurs to form 
deposit 


To 


Dimensionless Particle Size 


D. 


Empirical Parameter of Fluid Flow 


a 


* 




Location of the measurement point along a 
characteristic curve. 


s 0 
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Distance Parameter between the point of interest and 
the measurement point in characteristic curve 
coordinates. 


(s-Sn) 


Empirical exponent used in scaling deposit thickness 
to compensate for entrainment 


p 


Porosity of the deposit 


* 



Sediment Transport and Fluid Flow Physics 

[0021] In an embodiment of the present invention, the behavior of a fluid flow, 
and the characteristics of the deposit that flow forms, are described by the following 
equations. It will be understood to those skilled in the art that the functional reference 
(x,y) for parameters which are functions of x and y has been omitted for convenience 
in the following, but that omission is not limiting on the discussion which follows or 
the method of the present invention. 

[0022] Conservation of flow volume; 



10 



8hu dhv A 

H = U 



(i) 



dx dy 

The terms hu and hv represent flow volume flux per unit length in the x and y 
directions, respectively. Equation (1) means that the flow volume flux into any region 
within the fluid flow is equal to the flow volume flux out of that region 



15 



20 



[0023] Conservation of sediment: 

dx dy 



(2) 



The term v^d t ) is the effective settling velocity of grains of diameter d t . The effective 
settling velocity is related to the actual settling velocity of grains in still water, but in 
addition contains correction factors which, for example, may adjust for the fact that 
the water is not still. The term v s (d t )G therefore represents the rate at which the 
sediment volume of the z-th grain size is settling out of the flow. As a result, equation 
(2) indicates that for each grain size volume fraction Q, the sediment volume flowing 
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into a region is equal to the sediment volume flowing out of that region plus the 
sediment volume that settles out of the flow within the region. 

[0024J Deposit thickness: 



1-0 , 



(3) 



5 Equation (3) simply states that the thickness z of the deposit at any xy point may be 
determined by summing the settling rates for each grain size bin, multiplying by the 
time duration of the deposit, and dividing by the volume fraction of the deposit which 
is sediment (and not interstitial water), that is, dividing by one minus the porosity, <j> , 

of the deposit. 

10 [0025] Settling velocity for grains of sediment in still water: 

The effective settling velocity relationship used in one embodiment of the present 
invention is based on the empirical still water settling velocity relationship of 
Dietrich, "Settling Velocity of Natural Particles, " Water Resources Research, Vol. 
18, 1982, pp. 1615-1626. Dietrich defines a "dimensionless particle size" 



15 



D Jtf>s~ Pw)pwZ d \ 



(4) 



20 



and expresses the still water settling velocity w s as, 



gf*{p s -Pw) 



Dl 
5832 



1/3 



, ZL < 0.05 



gM(Ps~ 



r a 0 + a, log A* + a 2 (logA, f y 
Pw eXP | + a 3 (log/),) + a 4 (log D^f J 



-|l/3 



, A, 2t 0.05 



(5) 



where, a 0 =-3.76715 , a, =1.92944, a 2 =-0.09815, a 3 =-0.00575, and 
a 4 = 0.00056 . 

25 [0026] In one embodiment of the present invention, the effective settling velocity 
is defined as follows: 
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(D 0 ) a w s (D 0 ) d,*D t 



(6) 



where the multiplicative factor of grain diameter to the power a is applied to 
Dietrich's settling velocity. In addition, a cutoff diameter Do is specified below which 
the effective settling velocity is considered constant. Note that other alternative forms 
of the functional dependencies of effective settling velocity on grain diameter may be 
employed within the scope of this invention. 



[0027] lh the effective settling velocity definition of equation (6), the power a is 
10 an empirical parameter specified based on the presumed nature of the fluid flow that 
resulted in deposition of the deposit of interest. Rare large fluid flows of high 
velocity (flood conditions) account for most of the sedimentary deposits which 
become part of reservoir formations. Such flows tend to be suspended load- 
dominated, and empirical studies indicate that a value of a = 0 is appropriate. On the 
15 other hand, more common fluid flows of modest velocity tend to be bedload transport- 
dominated, and appear, from empirical studies, to have a settling velocity which has a 
much weaker dependence on grain diameter and a value of a = -1.5 appears most 
appropriate. 

20 [0028] Furthermore, empirical studies indicate that in deep water cases where the 
depositing flow persists a short time, grain size predictions are improved if a lower 
limit on the effective settling velocity is assumed. These studies suggest that a cutoff 
diameter, Do, of 25 microns be used. It is believed that a significant volume of these 
smaller grains are suspended in a cloud above the main body of the flow and do not 

25 follow the same flow lines as the larger sediment fractions. The need for an effective 
settling velocity that does not change with grain diameter below a threshold diameter 
suggests that these smallest grain size fractions form a well-mixed cloud of sediment, 
where the size distribution does not change laterally. In cases where the flow persists 
long enough that the sediment deposited by the flow is much greater in volume at all 
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ft 

points of interest than the sediment deposited by the suspended cloud, the cutoff 
diameter, Dq, should be 0 microns. 

[0029] Based on these principles, one embodiment of the method of the present 
5 invention can be derived as follows. First, using the rule for the derivative of a 
product on equation (2) and simplifying the result using equation (1) gives: 

hu^+hv^ = -v s (d t )C t (7) 

[0030] As will be understood to those skilled in the art, equation 7 is a linear 
partial differential equation of first order in Q. In the second step, this equation may 
10 be solved by the method of characteristics, as more fully described by Garabedian, 
Partial Differential Equations, 1964. This solution involves the introduction of a 
variable s, which may be defined as a measure of distance along a characteristic curve 
of the fluid flow, as follows: 

^- = hu, & = hv (8) 
ds ds 

15 Equation (8) allows equation (7) to be written as: 

^ (9) 

[0031] The physical significance of equation (9) is that the rate of change of grain 
volume fraction along a characteristic curve is proportional to the rate at which the 
sediment volume is settling out of the flow. This relationship holds to the extent that 
20 deposition rate may be considered independent of flow velocity. 

[0032] Note that the characteristic curves which result from introduction of the 
variable s will not necessarily be coincident with flowpaths of the fluid flow. 
Equation (8) involves the product of the height of the fluid flow, h, and scalar velocity 
components, u and v. As will be observed in Fig. 2, height h varies with location 
25 within the fluid flow, based for example on the height of the deposit z or flow velocity 
changes. These and other similar considerations are included in the flow volume flux 
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terals on the right hand side of equation (8) and cause the characteristic curves s to not 
necessarily be identical to fluid flowpaths. Note that the use of characteristic curves 
in the present invention are for calculational convenience, but are not a limitation of 
the present invention. 

5 [0033] Third, integrating equation (9) gives: 

C, W = C, (s 0 )exp[- v, {d,Xs - s 0 )] (10) 



[0034] Equation (10) means that the grain size volume fraction at a point s along a 
flowpath is related to the grain size volume fraction at an initial point s 0 multiplied by 
10 a decaying exponential term involving the effective settling velocity and the 
difference between s and s 0 . 

[0035] Fourth, combining this result with the deposit thickness relationship in 
equation (3) gives 

z M=r*lZ v.feJCifeM-v^Xf-^)] (11) 

15 [0036] Finally, it will be understood that, for a given location s, the grain size 
distribution G t {s) of grains in the deposit can be defined as: 

M-Satom d2) 



Description of Example Embodiments 

20 

[0037] This five step derivation allows a simple three step process, as depicted in 
Fig. 3, to be followed for predicting the grain size distribution at any location of 
interest in the deposit. First, the input parameters are specified, 31. This involves the 
determination or estimation of the thickness and grain size distribution at a single 
25 location in the deposit, for example from well data. For reference, this location, s 0> 
has a known deposit thickness z(s 0 ) and deposit grain size distribution G t {s 0 ). In 
addition, the deposit thickness z(s) at the location, s, of interest must be specified, or 
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estimated, for example from seismic data. Second, the distance parameter is 
computed, 33. The following equation is used to solve for the distance parameter (s - 
s 0 ): 

*W = z ( s o )£ G, (f o )exp[- v s (d t X* - «, )] (13) 

i 

5 Equation (13) results from combining equation (12) and the known values of z(s 0 ) and 
G t {s 0 ) with equation (11). Third, this result is used to determine G t {s) 9 35. In this 
example, G,(s) is determined from the following equation: 

G / W = 4r G /fc)exp[~v,(rf / X^ 5o )] (14) 
z[s) 

Equation 14 derives from substituting equation (10) into equation (12), and separately 
10 substituting the outcome of evaluating equation (12) at location s 0 into the resulting 
equation. 

[0038] This three step analysis process allows the grain size distribution G ( at 
every location s within the reservoir to be determined, solely on the basis of the three 
known or assumed input parameters: deposit thickness z(s Q ) 9 grain size distribution 
15 Gt{s 0 ) 9 and deposit thickness z(s). 

[0039] The assumptions embodied in equations 1, 2, and 3, along with the 
assumption that the effective settling velocity does not depend on the flow velocity, 
mean that each location within the deposit having a given thickness, referred to herein 

for convenience as Z* , will have a corresponding grain size distribution, referred to 

20 herein for convenience as G* . Furthermore, all locations within the deposit having 

the same thickness Z* will have the same grain size distribution G* . This thickness 
to grain size distribution relationship is the same along all characteristic curves (and 
hence over the entire area of the deposit), and is indexed by the single parameter 

(s-s 0 ). As a result, equations (13) and (14) involve calculations which are 

25 independent of whether the location of interest falls on the same characteristic curve 
as does the reference location . It should further be understood that the "distance 
parameter" (^-^ 0 ) is an effective distance in the characteristic curve coordinate 
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system and is not a physical distance between the reference location and the location 
of interest 

[0040] Noise in the measurement of grain size distribution and the physical 
approximations implicit in this method may limit the ability to extrapolate grain size 
5 distribution to locations where the thickness is several or more times larger or smaller 
than the measurement location. In such cases a modification of the method may be 
applied to achieve superior predictions. In this modification, grain diameters for at 
least one percentile of the grain size distribution (e.g., the "p50" grain diameter) are 
determined at the measurement location. The grain diameters for the same percentiles 

10 and an associated thickness are recalculated for a value of s-sq different from zero. 
The grain diameters corresponding to the same percentiles for any deposit thickness 
are then determined by a fitting a power law relationship to the grain diameters as a 
function of thickness at the measured location and the calculated location. For 
example, if the deposit thickness at the measurement location is z(so) and the grain 

15 diameter of the Mh percentile of grain size is ptfao) at the measurement location, and 
the projected thickness for a value of s different from sq is z(s), and the grain diameter 
of the Mh percentile of grain size is pri(s) 9 then the grain diameter at the Mh 
percentile of grain size at a different location s f with deposit thickness z(s r ) can be 
determined by fitting a power law relationship through the two grain diameter and 

20 thickness points. This is most easily accomplished by a linear extrapolation based on 
the logarithms of the diameters and thicknesses, according to equation 15, 



[0041] Turbidity currents in deep water entrain surrounding clear, stationary 
water. This entrainment modifies the relationship between thickness and grain size in 
25 the deposits of such currents. A simple modification to the inventive method 
compensates for the effect of entrainment, allowing the inventive method to be 
applied to deposits where entrainment is significant. To perform this modification, 
the measured deposit thickness z(s) at the designated location is first scaled to form 
z'(s) according to equation 16, 






WO 2004/083896 



PCT/US2004/005019 



-14- 



z\s) = z{s) 



£(£) 

A s o)_ 



(16) 



where z(?o) is the deposit thickness at the (second) location where the grain size 
distribution is measured. The exponent {$ has been determined empirically to be -0.3. 
5 The calculations of equations 13 and 14 are then applied as before, but using the 
scaled thickness z'(s) in place of the measured thickness z(s). Other choices of p are 
within the scope of this invention and may be used to partly compensate for other 
physical effects, such as the effect of flow velocity on the settling velocity of grains 

■ 

and the effect of porosity differences on measured thickness throughout the deposit. 
10 Other mathematical manipulations which achieve a similar adjustment to the 
predicted relationship between thickness and grain size distribution are also within the 
scope of this invention. 

[0042] The inventive method can be applied to composite sand bodies, which are 
composed of many smaller sand bodies. In this case, the inventive method operates 

15 on the total thickness of the composite body and the averaged grain size distribution 
within the composite body at the specified plan-view points. In one embodiment, the 
grain size averaging is performed by summing the grain size distributions within each 
smaller sand body at a specified plan-view point, weighted by the fractions of the total 
thickness which each smaller body represents. Other methods for computing an 

20 average grain size distribution for a composite sand body are also within the scope of 
this invention. 

[0043] In the inventive method, the predicted rate of change in grain size with 
thickness depends significantly on the sorting of the grain size distribution. Broader 
distributions generally result in predictions of more rapid changes in grain size with 
25 thickness. When the average grain size distribution of a composite sand body is 
calculated as the thickness-weighted average of the grain size distributions in the 
bodies which compose it, the average distribution tends to be broader than the 
distribution in most of the individual bodies. In this case, the inventive method would 
tend to predict too rapid a change of grain size with thickness, unless this effect is 
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corrected by adjusting the effective settling velocity by choosing a < 0, typically -1.5, 
in equation (6). In the alternative, the adjustment can be made by choosing a more 
negative value of P in equation (16). 

Typical Results: 

5 [0044] Figures 4 and 5 depict the results of a first application of an embodiment 
of the present invention. Figure 4 depicts results from a laboratory test in which 
sediment-carrying water was emitted at 16.5 cm/s with a 3.5% sediment concentration 
by volume from 2 inch diameter pipe 20 for five and a half minutes. Deposit 22 
formed in front of the pipe. The thickness of the deposit is indicated by the relative 

10 shading of the image. Distance measurements in the x and y directions are in units of 
centimeters. A grain size distribution sample was taken at location 24 and used 
together with the measured deposit thickness at each of the locations 25 indicated 
with an X to predict the grain size distribution at those locations. In addition, actual 
grain size distributions were also measured at the locations 25 indicated by X's for 

15 comparison with the prediction. Figure 5 shows the comparison of the measured to 
the predicted grain size distributions. Figure 5 plots deposit thickness on the 
horizontal axis in units of centimeters against the grain diameter on the vertical axis in 
units of microns. Three examples of predicted grain diameters are shown, 
corresponding to the tenth percentile (P 10) in the predicted grain size distribution, 26, 

20 the fiftieth percentile (P50), 28, and the ninetieth percentile (P90), 30, although it 
would be understood that the present invention predicts the entire grain size 
distribution G. These predictions are shown by the respective lines, with the 
measurements corresponding to each percentile shown by the corresponding symbols 
32, 34, and 36, for the respective grain size distribution percentiles. The symbols 32, 

25 34, and 36 with a grey-shaded border represent the measurement reference values 
used in the present method. In this example, the analysis assumptions included a = 0, 
Do = 25 microns, and p = 0. As can be observed, the measured values correspond 
well to the predictions resulting from an embodiment of the present invention. 

[0045] A second example of an application of an embodiment of the present 
30 invention is depicted in Figures 6 and 7. Figure 6 shows an aerial view of the Wax 
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Lake Delta in Southern Louisiana. X markings 50 on Figure 6 indicate locations 
where core measurements were taken. The vertically averaged grain size distribution 
was calculated at each core location, and the thickness of the delta was also measured 
at the same locations. Data from core 52 was used as the reference point, and the 

5 average grain size distributions were predicted at core locations 50. Figure 7 is 
similar in format to Figure 5, with the P10, P50, and P90 grain diameters for the 
predicted distributions shown as lines, 54, 56, and 58 respectively, and the measured 
values shown as symbols, 60, 62, and 64. Symbols 60, 62, and 64 with a grey-shaded 
border are the measurement reference values corresponding to core 52 in Fig. 6. In 

10 this application a = 1.5, Do = 0 microns, and p = 0 were assumed. As with Fig. 5, the 
measurements show good agreement with the predictions from the embodiment of the 
present invention. 
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CLAIMS 

What is claimed is: 

1 . A method for predicting the grain size distribution at a designated location 
within a water-lain sedimentary deposit, said method comprising: 

5 (a) determining the vertical thickness of said sedimentary deposit at said 

designated location; 

(b) determining the vertical thickness of said sedimentary deposit and the 
grain size distribution in said sedimentary deposit at a second location within said 
deposit; 

10 (c) determining a distance parameter corresponding to said designated 

location and said second location; and 

(d) using said distance parameter, said vertical thicknesses, and said grain 
size distribution at said second location to calculate said grain size distribution at said 
designated location. 

15 2. The method of claim 1 wherein said calculation involves an exponential 

dependence on an effective settling velocity of grains within a fluid flow from which 
said deposit was formed. 

3. The method of claim 2 wherein said effective settling velocity is a function of 
grain size. 

20 4. The method of claim 3 wherein said function of grain size involves a constant 
minimum diameter threshold. 

5. The method of claim 2 wherein said effective settling velocity involves an 
empirical relationship to the nature of said deposit-forming fluid flow. 

6. The method of claim 1 wherein said distance parameter involves use of a 
25 characteristic curve parameter. 
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7. The method of claim 6 wherein said characteristic curve parameter is defined 
in relation to a flow volume flux parameter corresponding to the fluid flow from 
which said deposit was formed. 

8. The method of claim 1 wherein said calculated grain size distribution is a 

5 function of the ratio of said thickness of said deposit at said second location to said 
thickness at said designated location. 

9. The method of claim 1 wherein said calculated grain size distribution is 
proportional to said grain size distribution at said second location. 

1 0. The method of claim 6 wherein said calculated grain size distribution is 
10 calculated from 

G/W=4r G '^>xp[-v,feX^^o)] 

z[s) 

wherein: 

Gi(s) is the fraction of total grain volume at location s in the deposit contained 
in the ith grain size bin 

9 

15 Gi(so) is the fraction of total grain volume at location so in the deposit 

contained in the ith grain size bin 

z(so) is the thickness of deposit at location so where the grain size distribution 
was also measured. 

z(s) is thickness of deposit at location s where the grain size distribution is to 
20 be determined. 

v s is the effective settling velocity of sediment grains as a function of grain 
diameter 
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dj is the nominal diameter of the grains in the ith grain size bin 

(s-so) is the distance parameter between the point of interest s and the second 
point so where the grain size measurement is available. 

1 1 . The method of claim 1 wherein said measured grain size distribution is the 
5 vertically averaged grain size distribution through a stratigraphic unit and said 

thickness at said designated location and said thickness at said second location are 
measured with respect to the top and bottom of the same stratigraphic unit, which may 
be comprised of one or more individual sedimentary deposits. 

12. The method of claim 1 wherein said calculated grain size distribution is 

10 compared with a measurement of grain size at said designated location to determine 
whether said designated location is within the same sedimentary deposit as said 
second location. 

13 The method of claim 1 1 wherein said calculated grain size distribution is 
compared with a measurement of grain size at said designated location to determine 
15 whether said designated location is within the same stratigraphic unit as said second 
location. 

14. The method of claim 1 applied with the designated location at every point 
within a sedimentary body such that the grain size distribution at every point within 
the sedimentary body is calculated. 

20 15. The method of claim 14 where the grain size distribution calculated for every 
point in the sedimentary body is weighted by the thickness at the respective point and 
averaged to determine the source grain size distribution for the system. 

16. A method for determining grain size from deposit thickness, where the 
relationship between grain size and thickness is extrapolated from at least one grain 
25 size point calculated according to method 1 . 
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17. The method of claim 16, wherein the extrapolation of grain size is based on a 
power law relationship between grain diameter at a given percentile of the grain size 
distribution and deposit thickness. 

18. The method of claim 1, wherein the deposit thickness at said designated location 
5 is first multiplied by a scale factor. 

19. The method of claim 18, wherein said scale factor is a power of the ratio of the 
deposit thickness at said designated location to the thickness at said second location. 
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